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ABSTRACT 



Context. Recent observational and theoretical works have suggested that the TuUy-Fisher Relation might be generalised 
to include dispersion-dominated systems by combining the rotation and dispersion velocity in the definition of the kine- 
matical indicator. Mergers and interactions have been pointed out as responsible of driving turbulent and disordered 
gas kinematics, which could generate TuUy-Fisher Relation outliers. 

Aims. We intend to investigate the gas kinematics of galaxies by using a simulated sample which includes both, gas 
disc-dominated and spheroid-dominated systems. We pay particular attention to the evolution of the scatter of the 
TuUy-Fisher Relation. We also aim at determining the gas-phase velocity indicator which better traces the potential 
well of the galaxy. 

Methods. Cosmological hydrodynamical simulations which include a multiphase model and physically-motivated 
Supernova feedback were performed in order to follow the evolution of galaxies as they are assembled. We analyse 
the gas kinematics of the surviving gas discs to estimate all velocity indicators. 

Results. Both the baryonic and stellar TuUy-Fisher relations for gas disc-dominated systems are tight while, as more 
dispersion-dominated systems are included, the scatter increases. We found a clear correlation between (r/Kot and mor- 
phology, with dispersion-dominated systems exhibiting the larger values (> 0.7). Mergers and interactions can affect the 
rotation curves directly or indirectly inducing a scatter in the TuUy-Fisher Relation larger than the simulated evolution 
since 2 ~ 3. Kinematical indicators which combine rotation velocity and dispersion velocity can reduce the scatter in 
the baryonic and the stellar mass-velocity relations. In particular, si.o = (l^ot + (J^)"'^ seems to be the best tracer of 
the circular velocity at larger radii. Our findings also show that the lowest scatter in both relations is obtained if the 
velocity indicators are measured at the maximum of the rotation curve. 

Conclusions. In agreement with previous works, we found that the gas kinematics of galaxies is significantly regulated 
by mergers and interactions, which play a key role in inducing gas accretion, outflows and starbursts. The joint action 
of these processes within a hierarchical ACDM Universe generates a mean simulated TuUy-Fisher Relation in good 
agreement with observations since 2 ~ 3 but with a scatter depending on morphology. The rotation velocity estimated 
at the maximum of the gas rotation curve is found to be the best proxy for the potential well regardless of morphology. 

Key words, galaxies: formation - galaxies: evolution - galaxies: structure 



1. Introduction 



The Tully-Fisher Relation (|Tullv &: FishedllQTTl hereafter 
TFR) is considered one of the most fundamental scaling 
relations of disc galaxies as it links two important proper- 
ties: their stellar (sTFR) or baryonic (bTFR) mass content 
and the depth of their potential well measured by the ro- 
tation velocity. Therefore, the study of the TFR at differ- 
ent redshifts (z) could help to trace the dynamics of the 
gas and the star formation activity within the dark matter 
haloes, and hence, to provide constraints for galaxy forma- 
tion models ( e.g. |Avila - Reese et al., .1998: Mo et al., 1998t 



lAvila-Reese et al.ll2008l) 
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The determinati on of the slope and zero point of the 
TFR at low (e.g. .Bell fc de Jondl200ll: iPizagno et al.l l2007l: 



iMever et all 126681: IGurovich et al.ll2010|) and intermediate 



and high redshifts (e.g. 'Conselice et al. 



2005 



2006; Atkinson ct al. 2007; Kassin et al. 



2007 



120081: ICresci et all i2009t IGnerucci et al.l 1201 iD has been 



Flores et al. 



Puech et al. 



the subject of numerous works. Although many authors 
suggest the existence of evolution, observational results 
have not converged yet, being also unclear if the evolu- 
tion is present in the zero point, the slope , or in both 
(|Vogt et al.lll996l 119971: iNakamura et al.ll2006D . Recent ob- 
servational results by iMiller et al.l ( 201l[) show evidence 
for a modest evolution of the sTFR by 0.04 ± 0.07 from 
z ~ 1 to z ~ 0.3 in the sense that, at a given rota- 
tion velocity, galaxies have lower stellar masses in the past. 
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According to the results of ICresci et"~all (|2009[ ). the sTFR 
seems to be already in place at z ^ 2 with a slope sim- 
ilar to the local one but displaced towards lower stellar 
masses by 0.41 ±0.11 dex. These trends are in general good 
agreemen t with predictions of hydrodynamical simula- 
tions (e.g. lPortinari k, Sommer-Larsenl2007l : lde Rossi et al.l 
l2010tl . At higher redshifts, it is not yet clear if the TFR 
is at place due to t he high scatter so far measured (e.g 
iGnerucci et "aIll201lD . The origin and evolution of the scat- 
ter of the TFR is also considered as an open problem. 
Observational results suggest that it might be driven by sys- 
tems which exhibit disturbed kinem atics as a consequence 
of ga l axy interactions and mergers pCannappan fc Barton 



2004 iFlores et al.l 120061: iPuech et al.l l2008t iKassin et al. 
20071 iCovington et al.ll20"l0r Ba sed on previou s studies of 
the TFR (jTuUv fc FouauelllQSSll . lWeiner et al.l (2006') pro- 
posed a new kinematical estimator generating a combined 
velocity scale, sk'- 
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where K is a, constant < 1 and a is the velocity disper- 
sion. The scale sk combines rotation (associated to order 
motion) and pres sure support ( assoc iated to turbulent and 
random motion). iKassin et al.l (|2007l ) reported that the ob- 
served sTFR resulting from the used of sq.s is capable of 
generating a unified relation for galaxies of all morpholog- 
ical types in their sample, including disturbed and merg- 
ing cases. This TFR based in sq.s was also found to be 
remarkably tight over 0.1 < z < 1.2. B y using hydrody- 
namical pre-prepared merger simulations, ICovington et al] 
(i201Q) showed that the kinematical indicator sq.s corre- 
lates with the total potential well of galaxies including 
baryons and dark matter. A close examination of the kine- 
matical scaling laws of galaxies led these authors to con- 
clude that the approp riate constant K should be ~ 0.5. 
ICovington et al.l (|2010f ) also found that the scatter of the 
TFR correlates with close encounters and mergers which 
suggests that the kinematics could be used to determ ine the 
merger stage of the galactic systems (see also Pedro sa et al.l 
|2008|) . Finally, these authors introduced the kinematical pa- 
rameter S = X So. 5 which seems to be a good tracer of 
the circular velocity in their simulations. Also, Weiner et al.l 
()2006[) proposed a modified version as si.o = ^/V^^^^+^. 

In the current paradigm for galaxy formation, systems 
assemble hierarchically in a bottom-up fashion so that more 
massive galaxies formed by the aggregation of smaller ones. 
In this scenario, mergers and interactions between galactic 
systems play a major role in the evolution of baryonic mat- 
ter within the dark matter potential wells, partially reg- 
ulating the gas infall and star formation process, which 
also influence the efficiency of feedback mechanisms. In this 
context, cosmological N-body/hydrodynamical simulations 
constitute an optimum tool to follow the evolution of the 
astrophysical properties of galaxies and may help to the 
understanding of the fundamental relations between them. 
Our simulations have been run with a version of GADGET- 
3 which i ncludes a physically-mo tivated Supernova (SN) 
feedback (jScannapieco et aLll2006D which is able to trigger 
mass-loaded gas outflows. This scheme does not require ad 
hoc mass dependent parameters but naturally regulates the 
strength of the SN feedb ack in haloes of different masses 
(jScannapieco et 311120081 ). 



Fig. 1. Distributions of disc-to-total mass ratios (Fdisc) 
estimated by using the gas (solid line) , stellar (dashed line) 
and baryonic (shaded area) mass in simulated galaxies at 
z — 0. All simulated galaxies have a surviving gaseous disc 
even in systems where stars are all forming an spheroidal 
component. 

In Ide Rossi et"aI1 (|2010D . we studied the sTFR and 
bTFR in simulated galaxies with very well defined- 
dominating gaseous discs (disc-to-total gas mass ratios 
D/T > 0.75) formed in a A-CDM universe. We found 
that both the sTFR and the bTFR are well reproduced 
by our simulations. Our SN feedback model works in a 
self-regulated way so that small systems are more largely 
affected than large ones. Our simul ated sTFR i s in ag ree- 
ment to that recently estimated bv iReves et al.l (|201l[ ). in 
the same mass range. While the sTFR is detected to have 
a bend in the slope for low stellar mass systems as a re- 
sult of the SN feedback, the bTFR does not show the same 
feature except at very high redshifts. So, in this previous 
work, we found that the bTFR does not store information 
on the impact of SN feedback for reasonable SN energies, 
except at very early times (z ~ 3). 

In this paper, we extend the work of de Ros si et al.l 
(|2010D studying the gas kinematics of all galaxies with 
surviving gaseous discs, even those systems dominated by 
spheroidal components, using the same set of cosmological 
hydrodynamical simulations. Particular attention is paid to 
the study of the origin and evolution of the scatter in the 
TFR considering the impact of mergers and interactions 
as well as the importance of gas inflows and outflows. The 
rotation curves are derived from the analysis of the sur- 
viving gaseous discs with the aim at determining the best 
gas-phase kinematical tracer for the gravitational poten- 
tial. We assumed that, even if the surviving gaseous disc is 
tenuous, it will trace the potential well in which it inhabits. 

In section |21 we introduce the numerical simulations 
studied in this work, the galaxy catalogue and the methods 
developed to perform the analysis of gas kinematics. In sec- 
tion|21 we present a discussion of the local TFR and various 
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Fig. 2. Gas fraction as a function of stellar mass for sim- 
ulated galaxies. Different symbols depict systems with dif- 
ferent values of F^^^^, as indicated in the figure. Symbols 
enclosed by squares correspond to galaxies with F^^^^ > 0.5. 
The black dot ted curve are t he an alytical fits to observa- 
tions given bv lStewart et al.l (|2009l ) at z = 0. 



kinematic gas indicators for galaxies of different morpholo- 
gies. In sectional we investigate the fingerprints of the hi- 
erarchical aggregation of the structure on the TFR-plane. 
The conclusions of the work are summarised in section [5] 



2. Numerical Experiments 

The simulation analysed in this work was performed by us- 
ing the chemical code GADGET-3, an update of GADGET-2 
optimised for massiv e parallel simulations of highly inho- 
mogcncous systems ijSpringel fc HernquistI l2003t ISpringell 
L2005). Our version of the code includes treatments for 
metal-dependent radiative cooling, stochastic star forma- 
tion, a multiphase model fo r the interstellar medium (ISM ) 
and a SN feedback scheme (jScannapieco et ani2005l [20061) . 
The multiphase model allows the coexistence of diffuse and 
dense gas phases, impro ving the description of the ISM 
(jScannapieco et al]|2006[ ). The SN feedback model is able to 
trigger mass-loaded galactic outflows without introducing 
mass-dependent parameters. 

The chemical e volution prescri p tion u sed in this code 

(|2001|) and, later on, 
for GADGET-2. The 
model describes the enrichment by Type II (SNII) and 
Type la (SNIa) Supern ovae according to the c hemi- 
cal yield prescrip t ions of IWooslev fc Weaver! (|l995l) and 
[Thielemann et al.l (|1993D . respectively. We adopted a stan- 
dard Salpeter Initial Mass Function with a lower and upper 
mass cut-offs of 0.1 Mq and 40 Mq, respectively. It is as- 
sumed that each SN event generates 0.7 x 10^^ erg. For 
SNIa, we assumed a time-delay for the ejection of mate- 
rial randomly chosen within [0.1, 1] Gyr. We assumed that 



was develo ped bv iMosconi et al. 



adapted bv IScannapieco et al. 



SNII evolved within the same time interval in which they 
are formed. Metals are distributed within the neighbouring 
ga s particles weig h ted by the smoothing kernel as proposed 
bv IMosconi et al.l (|200l() . 

We use numerical hydrodynamical simulations consis- 
tent with a A-CDM universe with n = 0.3, A = 0.7, ilb = 
0.04, a normalisation of the power spectrum of erg = 0.9 
and Ho = 100 h km s-^Mpc"^ with h = 0.7. Although this 
set of parameters differs slightly from WMAP-7, the vari- 
ations have no important consequences on our results. We 
simulated a typical field region of the Universe in a comov- 
ing cubic volume of 10 Mpc h^^ side length. The simulation 
studied in this paper was run with 2 x 230'^ particles obtain- 
ing an initial mass of 9.1 x 10^ MqH'^ for gas particles, and 
a mass of 5.9 x 10^ M(T,h ~^ for the dark matt er particles. 
It is t he so-called S2 3 in [p j Rossi et al.l ()2010f ). 

In lde Rossi et al.l (j2010[ ). we showed that the dynamical 
properties of the galaxies in this sample are robust against 
numerical resolution and small changes in the parameters 
of the SN-feedback model. We also found that this sim- 
ulation predicts correlations between the dynamical prop- 
erties of galaxies in general good agreement with observa- 
tions, providing a suitable sample to study the origin and 
evolution of the TFR in a cosmological context. In par- 
ticular, we found that the stellar mass fractions of galax- 
ies residing within dark matter haloes are in good agree- 
ment with estimations derived from semi-empi rical mod- 
cls [3 bas ed on the observ ed l uminosity function (iGuo et al.l 
,2010; M oster et aLll2010[ ). In Ide Rossi et al.l (|2010f ). it was 
also proven that the adopted SN feedback model is able 
to reproduce successf ully the observed bend of the sTFR 
at ~ lOOkms"^ (e.g. iMcGaugh et al.l[2000l: lAmon'n et cH 
120091 : iTorres-Flores et al.ll20lH) as a consequence of its high 
efficiency at regulating the star formation p rocess of galax- 
ies in haloes with shallow potential wells (jde Rossi et al.l 
I2010D . On the contrary, the bTFR shows no SN feedback 
imprinted feature at low redshift. 



2.1. The galaxy catalogue 

The virialized structures in the simulated box are identified 
by a standard friends-of-friends technique and the substruc- 
tures within each dark matter halo a re then individualise d 
by using the SUBFIND algorithm ([Springel et al.l |2001D . 
In order to diminish numerical problems, we only analysed 
substructures with a total number of particles iVgub > 2000 
from z ~ 4. This gives a total number of galaxies which 
ranges from 150 at z = 3 to 309 at z = 0. 

The dynamical properties of the simulated galaxies have 
been estimated within the baryonic radius i?bar, defined as 
the one which encloses 83 per cent of the baryonic mass of 
the system. In these simulations, a typical Milky- Way type 
galaxy is resolved with around 10^ total particles within 
i?bar- It is also worth noting that for our simulated galaxies 
more than 90% of the gas inside i?bar have T <Tc — 8x 10^ 
K, with Tc being the critical temperature used by th e 
multi-phase model of the ISM (jScannaoieco et aLll2006t ). 



^ Note that those semi-empirical studi es used different Initial 
Mass Functions (Chabrier/Kroupa) to Ide Rossi et all l|2010l . 
Salpeter) and, hence, for the same luminosity implies some 40% 
less mass locked in stars. Even taking into account these dif- 
ferences, the simulated stellar mass fractions remain within the 
observational predicted range. 
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Fig. 3. Upper panels: Rotation curves for two typical simulated galaxies Gl and G2 (left and right panels, respectively) 
at z = with different fractions of gas supported by rotation (-Ffig^.)- The total circular velocity is shown with solid 
line. The circular velocities contributed by the gas, stellar, baryonic and dark components are represented with dotted, 
dot-dashed, dashed and long-dashed lines, respectively. The mean rotation curve (diamonds) and the mean dispersion 
velocity curve (squares) for the gas-phase are also plotted. The vertical lines indicate the baryonic radius (i?bar) and 
the radius where the mean rotation velocity reaches a maximum value (i?max)- Lower panels: Gas (solid lines) and 
stellar (dashed lines) mass in the spheroidal (red) and disc (blue) components as a function of radius, normalised to the 
corresponding total gas or stellar mass within S-Rbar- The vertical line depicts i?bar- 



Therefore, for the sake of simplicity, we did not distinguish 
between a cold and hot component in the analysis of gas 
dynamics as the mean behaviour is almost entirely deter- 
mined by the cold gas. 

2.2. Kinematic analysis of simulated galaxies 

For a given simulated galaxy, we calculated its circular ve- 
locity as T4irc = \/GM{r)/r, where M{r) is the total mass 
enclosed within the radius, r. In equilibrium, for baryons 
dominated by rotation, the rotation velocity V^-ot ~ K;irc, 
which is a measure of the potential well of the system. For 
each gas (stellar) particle i within a given galaxy, we defined 
the rotation velocity (l^ot,i) as its tangential velocity on the 
plane perpendicular to the total gas-phase (stellar-phase) 
angular momentum Jg (Jg) of the system. We estimated 
the averaged values of the velocity components for particles 
at equally-spaced radial bins in the gas and in the stellar 
components. Then, we calculated the dispersion velocity 



(aj) for each parti c le i in side these bins. Finally, following 
IScannapieco et all (|2005D . we used the ratio between the 
dispersion and the rotation velocity of each baryonic par- 
ticle {(Ti/Vj-ot^i) to define the disc and spheroid component 
of the given galaxy. If ai/Vrot,i < 1, we considered that the 
particle i belongs to the disc component, otherwise it was 
assumed that it belongs to the spheroid component Q. Note 
that our definition implies that the spheroid component is 
made of all particles dominated by dispersion so that it 
includes both the bulge and the inner stellar halo. 

In order to characterise the morphology of simulated 
galaxies, we estimated F^^^^ (^disc) ^he mass fraction of 
the total gas (stellar) mass within 3i?bar which is supported 
by rotation (i.e. the D/T ratio for the gas or stellar phase, 
respectively). We found that, in general, our sample of sim- 



^ We checked that this criterion produces simil ar results to 
those obt ained by applying the rn eth ods used by lAbadi et al.l 
(|2003h or IScannapieco et all (|2009l 'l or lTissera et al.1 (|2012D . 
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Fig. 4. Simulated sTFR (upper panels) and bTFR (lower panels) at z = by using Kot estimated at different radii as 
indicated in the labels. The green circles depict systems with F^^^^ < 0.75, while the blue squares represent systems with 
-^disc — O-^S. The black solid lines indicate the mean relations estimated by using Vcirc measured at the corresponding 
radii and their standard deviations. 



ulated galaxies is constituted by systems with dominating 
stellar spheroids and thick stellar discs while the gas mass 
fractions associated to the spheroid or the disc are very as- 
sorted. This behaviour can be appreciated more clearly in 
Fig. [1] where we show the histograms of F^-^^^ (solid line), 
-^disc (dashed line) and a similarly defined fraction for the 
baryonic mass (-Fdisd shaded area) for simulated galaxies 
at = [3. We can see that the percentage of the stellar 
mass dominated by rotation does not exceed 60% for any 
of the analysed galaxies while, in the case of the gas com- 
ponent, we obtained percentages ranging 20% to ~ 100%. 
When considering the baryonic phase as a whole, the per- 
centage of mass dominated by rotation exhibits an inter- 
mediate behaviour with percentages ranging between 10% 
and 80%. Most of the simulated galaxies in our sample have 
an important spheroidal component and all of them have a 
surviving gaseous discs from which it is possible to estimate 
rotation curves. Our goal here is to analyse if, regardless of 
galaxy morphology, the gaseous discs can defined a TFR 
consistent with observations and if the corresponding ro- 
tation curves are able to trace the potential wells of their 
galaxies. 

As shown in Fig. [2l the /gas of our simulated galax- 
ies decreases with stellar mass, consistently with observa- 
tions, although, at a given stellar m a,ss, they have lower 
gas fractions than observed ones (e.g. IStewart et al.ll2009t 
iReves et al.ll201lD as a consequence o f the very efficient SN 
feedback assumed in this model (e.g. lde Rossi et aLll2010D . 
We also appreciate that simulated galaxies with dominant 
baryonic discs {F^^^^ > 0.5) are biased to lower stellar 
masses and high gas fractions indicating that the baryonic 
discs are formed by significant gas masses. In this figure, we 
also highlighted the simulated galaxies according to the pa- 

^ ^diac' ^disc ''■nd -Fdisc equivalent to the commonly used 
D/T ratio for the gas, stellar and baryonic components. 



rameter (F|,^ < 0.5, 0.5 < F^,, < 0.75, F^,^^ > 0.5) 
to show the large variety of gas disc-type galaxies coming 
out of the simulations. While there are not massive stellar 
galaxies (M, ^ lO^'^MQh"^) with important stellar disc 
components, the left-over gas mass is mainly dominated by 
rotational motions. We have no system dominated by stellar 
discs as indicated by the lack of population with F^^^^ < 0.5 
and i^b,, > 0.5. 

Normal spiral galaxies have their gaseous and stellar 
disc components aligned. Cases of misalignment are gener- 
ally considered the product of interactions or rece nt mergers 
which perturbed the kinematics of the systems (jSnaith al.l 
I2OI2D . Therefore, in order to have a more clear picture of 
the kinematics of baryons in our sample, we calculated the 
cosines of the angles (a) determined by the angular mo- 
menta of the gas discs ( Jg) and those of the stellar ones 
(.7*). In the case of galaxies with dominant baryonic discs 
(-^disc — 0-5), we estimated these cosines taken into account 
only the mass associated to the disc component while, for 
systems with dominant baryonic spheroids {F^^^^ < 0.5), 
we considered the whole baryonic mass. Our results show 
that cos(a) > 0.95 for all systems with F^^^^ > 0.5, showing 
that gaseous and stellar discs are very well aligned, consis- 
tently with disc-like galaxies. In the case of galaxies with 
-^disc ^ '^•5' "^g J* are not always aligned, forming an- 
gles larger than 7r/4 in some cases, which reveals different 
levels of disturbance in the baryonic kinematics. For fu- 
ture discussion, we will define three subsamples: SI, consti- 
tuted by systems with dominant baryonic discs, which are 
more similar to real spiral galaxies; S2, formed by systems 
with dominant baryonic spheroids and aligned gas discs 
(^disc < 0-5 and cos(a) > 0.7), and S3, including galax- 
ies also with dominant baryonic spheroids but misaligned 
gas discs (iodise < ^-^ ^^d cos(a) < 0.7). The percentages 



5 



De Rossi et al.: Gas kinematics of galaxies 



2.2.1. Rotation curves 
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Fig. 5. Kot/V^irc and cr/Vrot as a function of the fraction of 
the gas mass supported by rotation and at different radii in 
the simulations at z = 0: i?max (dot-dashed line), 0.5i?bar 
(solid hue), i?bar (dotted line) and 1.5i?bar (dashed line). 
Error bars correspond to the standard deviations. It is in- 
teresting to note that V^ot at i?max can approximate V^irc 
with good accuracy (> 90%) for all simulated systems. 



of galaxies corresponding to SI, S2 and S3 are 33%, 31% 
and 36%, respectively. 

As the main goal of this work is the analysis of the 
gas kinematics in the disc component, we will characterise 
simulated systems by using only -Ffjg^ but, when necessary, 
we will also refer to the subsamples SI, S2 and S3 in order 
to get a more complete picture. Note that in lde Rossi et al.l 

> 



we restricted our study to galaxies with F^^ 



disc 



0.75 so that the gaseous disc components were important 
with well-behaved rotation curves. As mentioned before, in 
this work, we extend that analysis to the whole sample of 
simulated galaxies with a surviving gas disc. 



In the upper panels of Fig. |31 we show the rotation and 
circular velocities for two typical simulated galaxies at z = 
0. The left panels exhibit the results for a galaxy with a 
stellar mass of A/, = IO^-^Mq/i-i and Ff;^^ = 0.9 (Gl), 
while in the right panels we display those of a galaxy with 
a stellar mass of Af, = IO^^MqH-^ and F^-^^ = 0.6 (G2). 
Gl belongs to the subsample SI (i.e. systems more similar 
to real spiral galaxies) and G2 to the subsample S2 (i.e. 
systems with dominant baryonic spheroids, with aligned gas 
discs). As can be seen, Gl has a more important gas-phase 
disc component with a rotation curve tracing the circular 
velocity remarkably well until i?bar, while in the case of 
G2, it drops around -Rmax, departing significantly from the 
circular velocity. Note that, in both cases, the decrease in 
the rotation curve is accompanied by an increase of the gas 
mean dispersion velocity, as expected. 

In the lower panels of Fig. [31 we can appreciate the gas 
(solid lines) and stellar (dashed lines) mass fractions associ- 
ated to the spheroidal (red) and disc (blue) components as 
a function of radius for Gl and G2. The gas (stellar) mass 
in each component is normalised to the total gas (stellar) 
mass enclosed by 3i?bar- We see that, in the case of Gl, the 
stars are more concentrated than in G2. 



3. The Tully-Fisher relation 

The kinematical analysis performed for Gl and G2 in sec- 
tion 12.21 was extended to our whole sample of simulated 
galactic systems at z = 0. We constructed the sTFR and 
bTFR by using the rotation velocity and the circular ve- 
locity estimated at different radii: i?max, 0.5i?bar, -Rbar and 
1.5i?bar, as shown in Fig. S] Ou r simulated TFR a grees 
very well with that reported by iReves et al.l ()201lD over 
the whole observed mass range (« 10^ — 10^^MQh~^). 
Our SN feedback model is able to reproduce the observed 
sTFR and bTFR even for these large stellar masses with- 
out the need to reso rt to other physical processes (see 
'McCart hv et al.l l2012l for different results). As shown in 
ide Rossi et al.l (|2010D . this SN feedback model can also pre- 
dict the observed bend of the sTFR at smaller stellar masses 
A'U < lO^Afo/i-i (e.g. lMcGaugh et al.|[2000l: lAmon'n et"all 
2009). Indeed, the good agreement of our simulated TFRs 
with observations of large and small stellar masses is the re- 
sult of the action of the self-regulated SN feedback adopted 
in our models, which is able to regulate the star formation 
process in systems of different potential wells without the 
need to introduce ad hoc mass-dependent parameters. 

In the case of the gaseous disc-dominated systems (blue 
squares in Fig. H > 0.75, 23% of the total sam- 

ple), the good agreement between the TFR obtained by 
using Vj-ot and that estimated by using Vdrc at different 
radii, shows that the former is a good measure of the 
total potential well. For these systems, we obtained that 
the scatter in the sTFR and bTFR decreases when Vrot 
is esti mated at ^ 0.5i?bari in agreement with observa- 
tions (jYeeorova fc Saluccill2007f) . Also, for these systems, 
we found that T^ot (-Rmax) ^ V r ot(-Rb ai) Consistently with 
the findings of iPersic" fc Saluccil (|1995ll . 

In the case of systems with more important gas 
spheroidal components (green circles in Fig. |4l F^^^^ < 
0.75), their rotation velocities tend to be smaller than the 
Vch-c at a given mass. This is accompanied by the increase 
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Fig. 6. Vi/Kirc at 2 = as a function of radius in the simulations for Vi: a (blue triangles), Vfot (red squares), so.5 (black 
crosses), si.o (green asterisks) and S (pink diamonds). Each symbol represents the averaged value over the entire sample 
of galaxies (upper left panel) and the subsamples SI, S2 and S3 defined in section [22] (upper right panel, lower left panel 
and lower right panel, respectively) . The corresponding standard deviations are also shown. Note that the radii are given 
normalised to i?bar in the main panels while, in the insets, they are given normalised to i?max- The horizontal solid line 
indicates where Vi agrees with Vcirc, while the dashed horizontal lines denote the range where Vi does not depart more 
than 20% from T4irc- 



in gas velocity dispersion. To quantify the behaviour of 
the kinematical estimators, in Fig. [S] we show Vrot/Ki:irc 
and a /V-cot as a function of F^^^^, measured at different 
radii (cr is estimated by considering all gas particles within 
a given radius). When the disc dominates the gas-phase 
component, V-cot can approximate Vcirc with good accu- 
racy (larger than ~ 90%). However, as decreases, 
K-ot/Kirc becomes smaller and cr/Vrot gets larger as the 
gas component is more dominated by dispersion, showing 



more important variations when measured at larger radii. 
For galaxies with gas component dominated by dispersion 
(-^disc < "^e found a mean value of cr/Vrot ^ 0.7 mea- 
sured at i?Tna,y, in agreement w ith recent observational re- 
sults bv lCatinella et all (|2012| ). 

We recalculated these trends for the subsamples SI, S2 
and S3 (see section [2?2|) . obtaining similar global results. 
In the case of cr/K-ot , regardless of the alignment between 
stellar and gaseous angular momenta, galaxies with dom- 
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inating baryonic spheroids have also more disturbed gas 
dynamics, implying higher ratios than those for systems 
with baryonic dominating disc, at a given F^^^^. 

As we note from Fig. [5l a very interesting result from 
our model is that Vrot at i?max seems to be the best proxy 
for the circular velocity as it shows no significant depen- 
dence on T4ot is always within 90% of Vdrc for all the 
systems. This trend is even better for gas systems domi- 
nated by rotation. We checked these trends for SI, S2 and 
S3, finding similar results. Regardless of the relative angu- 
lar momentum orientation between stellar and gas compo- 
nents, V^ot at -Rmax provides, on average, the best indicator 
for the total circular velocity, although with high standard 
deviations for systems with F^^^^ ^0.5. 



As d iscussed in the Int roduction, iKassin et al.l (|2007[ ) 
and ICovi ngton et al.l (|201Cl[ ) proposed new kinematical es- 
timators based on observational and numerical results, re- 
spectively, to account for the disordered motions in the gas 
component with the aim at explaining the origin of the 
scatter in the TFR. We tested these estimators in our cos- 
mological simulations. In Fig. [6l we can appreciate Vi/Vdrc 
at 2: = as a function of radius where Vi could be either a, 
Kot, So. 5, si.o or iS*. Each symbol represents the averaged 
value over the entire sample of galaxies (upper left panel) 
and the subsamples SI, S2 and S3 defined in section [521 
(upper right panel, lower left panel and lower right panel, 
respectively). The corresponding standard deviations are 
also shown. Note that the radii are given normalised to i?bar 
in the main panels while, in the insets, they are given nor- 
malised to Rmax- The horizontal solid line indicates where 
Vi agrees with Vdrc, while the dashed horizontal lines de- 
note the range where V does not depart more than 20% 
from Vdrc- By analysing the main panels of this figure it is 
clear that a is lower than Vdrc at almost all analysed radii 
with the only exception of the inner parts of the galaxies. 
At radii larger than 0.5i?bar, f/V^irc remains almost con- 
stant at « 0.55 if we consider the whole sample of galaxies, 
« 0.50 in the case of SI and « 0.60 for S2 and S3. These 
findings are consistent with the fact that SI is comprised 
by more disc-like galaxies and hence, the dispersion veloc- 
ities tend to be smaller. For the whole galaxy sample, our 
results indicate that K-ot can be considered a proxy for Vdrc 
in the range 0.2i?bar < r < 0.5i?bar- Similar trends can be 
appreciated for the subsamples S2 and S3 but, in the case 
of SI, this range can be extended towards 0.7i?bar, approxi- 
mately. By comparing V^irc with so.5 and S, we see that the 
former underestimates V^irc, while the latter overproduces 
it for the whole sample, and for SI, S2 and S3. Nevertheless, 
in the case of SI, S constitutes a better representation of 
Vcirc as a consequence of the lower dispersion velocities in 
this subsample. The analysis of the rotation and dispersion 
velocity curves of our simulated galaxies suggests that the 
best proxy for V^irc at z = is si.o, which approximates 
remarkably well the simulated Vdrc at r > 0.5i?bar for our 
whole sample and each one of the subsamples. 

When analysing the behaviour of Vi/Vdrc as a function 
of R/Rmax (insets of Fig. (5]), we note that for the whole 
sample of galaxies the scatter in the vertical axis increases 
significantly, mainly at R/Rmax 1- We can see that the 
main contribution to this scatter comes from galaxies in S2 
and S3, while in the case of SI the scatter is quite small. For 
S3, the scatter is even larger than for S2 as these systems 
exhibit a more disturbed baryonic kinematics (misaligned 
gas discs). It is interesting to note that, in these simulations. 
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Fig. 7. Ratio between the standard deviations of V and 
V^irc (AVi/AVcirc) at z = for different bins of stellar 
(M,) and baryonic (Mb) mass. Each bin is constituted by 
35 galaxies of similar stellar or baryonic masses. Results 
are shown for V: Vrot (upper panels), so.5 (middle panels) 
and Si.o (lower panels). The kinematical indicators are es- 
timated by evaluating the gas kinematics at different radii 
in the simulations: i?max (dot-dashed line), 0.5i?bar (solid 
line), i?bar (dotted line) and 1.5i?bar (dashed line). 



the rotation curves of spheroid-dominated systems tend to 
be more affected in the inner parts of the galaxy (R ~ 
^max) than outwards (i? ~ -Rbar)- It is also worth noting 
that for the three subsamples, on average, Vrot ~ Vdrc at 

R ^ Rmax- 

In Fig. [71 we plotted the ratio between the standard 
deviation of the TFR by using Vrot, so.5 and si.o and the 
standard deviation by using Vdrc- The ratio is analysed 
as a function of stellar and baryonic mass. We can ap- 
preciate that by combining Vrot and a in the definition of 
the kinematical indicator sq.s or si.o, the dispersion in the 
sTFR and bTFR at a given mass is considerably reduced. 
By analysing these trends for SI, S2 and S3 separately, we 
found that this scatter is mostly generated by galaxies in S2 
and S3 (but mainly in S3) given the lower mass fractions as- 
sociated to their surviving discs. Nevertheless, all kinemat- 
ical indicators including Vrot leads to the tightest relation 
if evaluated at i?niax- Hence, V-ot at i?rnax not only seems 
to be a very good representation of the potential well of 
simulated systems, regardless of the morphology (Fig.[Slini) 
but also, the scatter of the sTFR and bTFR tends to be re- 
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duced if the velocity indicators are evaluated at this radius. 
This fact may be suggesting that -Rmax can be considered a 
characteristic radius for the gas kinematics of galaxies, at 
least in the case of the potential wells reproduced by these 
simulations. From Fig. [3 we also appreciate that both, so.5 
and Si.o reduces the scatter of the TFR by similar amounts 
but, in the literature, most authors use S0.5. In the case 
of these simulations, we found si.o to be a better tracer 
of the gas kinematics as it is also a good proxy for Vch-c 
(Fig. El). Finally, by comparing the relative deviations of 
the sTRF and bTFR (left and right-hand panels of Fig.©, 
it is evident that the bTFR is tighter than the sTFR for the 
analysed mass range everywhere, but at the lowest masses. 
The less massive systems are also the most affected by nu- 
merical resolution so the trends at the low-mass end should 
be confirmed with higher numerical resolution simulations. 

4. The effects of galaxy assembly in the TFR-plane 

In order to understand the evolution of the gas kinematics 
of galaxies and its implications for the origin of the scat- 
ter in the sTFR and bTFR, we analyse the evolution of 
galaxies on the TFR-plane as a function of time and in re- 
lation to important events such as mergers. To accomplish 
this purpose, we reconstructed the assembly history of the 
complete sample of simulated galaxies from their forma- 
tion time to z = and explored which events disturb their 
rotation curves leading to outliers on the TFR-plane. We 
designed our own algorithm to follow the main progenitor 
branch for each simulated galaxy. At each available time, 
we define the main progenitor as the most massive baryonic 
substructure. All the other smaller systems in the tree are 
considered satellites which will eventually merge onto the 
main branch. 

In order to evaluate the effects of the galaxy-galaxy 
interactions, for a given galaxy, we determined the total 
mass within i?bar (A^dyn") ^^'^ the mass of all its neigh- 
bours {M^°'^^^) located at a distance smaller than a cer- 
tain radius which can be taken either as the virial radius 
(i?vir) or as twice its value. We quantified the strength of 
interactions by the ratio PP^^^ = E ^^dyjf''^A^dyn''- This 
ratio is measured at both distances in order to have an 
idea of the satellite distribution with respect to the main 
progenitor (for example if most of them are within i?vir 
or further away). To quantify the strength of mergers, the 

ratio i^n^ergor = E ^dy^n Wyn" defined, whcrC M^^f^ is 

the mass of the satellites which were accreted between two 
available times. 

For the main progenitor, we estimated the gas fraction, 
its baryonic mass normalised to the value at z = and 
-^disc ^ function of lookback time in order to quantify 
the building up of the galaxy as well as the role of outfiows, 
as we will discuss in next section. We also calculate the star 
formation rate efficiency (SFRe) in the progenitor galaxy as 
the ratio between the star formation rate and the gas mass. 
We will use the ratios K-ot/K:irc and si.o/Vi^irc to quantify 
the disturbance of the rotation curve along the evolutionary 
paths. All these quantities will be estimated at i?bar- 

4.1. Formation histories of typical galaxies 

We analysed the evolution of the quantities defined above 
for the main progenitor branch of the local simulated galax- 



ies as a function of lookback time. In order to illustrate the 
relation between interactions and mergers and disturbances 
of the rotation curve, we selected six typical galaxies out of 
a sample of 309 galaxies at z = 0, as shown in Fig. [51 

The galaxy in the upper-left panel of Fig.|Slis one of the 
more massive galaxies at z = exhibiting a lower scatter in 
the local simulated TFR. Note that, as expected, Vi-ot/V^irc 
is approximately 1 at z ^ while, si.o/Kirc tends to be a 
little larger by ^ 0.1 dex. By analysing its formation his- 
tory, we see that the main progenitor has not experienced 
important interactions or merger events during the last 5 
Gyr. Its SFRe has also decreased by more than an order 
of magnitude in the same time interval as a consequence of 
the decrease of the gas fraction. Therefore, the galaxy at 
z = has reached a passive stage of evolution and stability 
which is consistent with its well-defined rotation curve and 
low scatter in the TFR. Note also that, although the gas 
fraction is small (~ 0.07) at z = 0, the gas mass within 
i?bar is « 3 X IO^Mq/i"^ so that it is possible to have a well 
resolved disc component. By exploring the main progeni- 
tor branch at early times, we can see that for a lookback 
time greater than 8 Gyr the main progenitor exhibited a 
disturbed kinematics. At this epoch, the ratio Kot/K:irc 
shows variations by ±0.1 dex over time intervals of around 
1 Gyr. From the inspection of 

-Winter and -^merger ^t that 

time, we can infer that this behaviour can be associated 
with interaction and merger events which lead to a high 
SFRe and exhausted the gas reservoir of the system. It is 
interesting that merger and interactio ns can generate eithe r 
an increase or decrease of K-ot/Vcirc (jPedrosa et al.ll2008h . 
With respect to si.o/Kirc, we obtained similar trends but 
with smaller and more smooth variations. 

The galaxy in the upper right panel is an example of 
an intermediate-mass galaxy with low scatter in the local 
TFR. This galaxy also exhibits a passive evolution within 
the last few Gyrs. Similarly to the previous example, we can 
see that at early times (^ 9 Gyrs ago), its main progeni- 
tor experienced important interactions followed by merger 
events {F^J^^f ^ 0.4, F^^™ ~ 0.3, F„,erger 0.25) which 
could be finked to a decrease in Vrot/K;irc ratio by around 
~ 0.5 dex over the corresponding period. Note also that 
^disc exhibits a minimum at that epoch. After these events, 
its disc component was reconstructed and the gas rotation 
velocity approached the circular velocity again. It is worth 
noting that si.o/K:irc ~ 1 during the whole formation his- 
tory of this galaxy, even during the epoch of mergers an 
interactions. As previously reported in the literature, these 
findings suggest that the combination of V^ot and cr in the 
definition of the kinematical indicator generates a velocity 
scale more robust against disturbances of the gas kinemat- 
ics. 

In the middle panels of Fig. [51 we compare the for- 
mation histories of two smaller galaxies with similar stel- 
lar masses of « IO^-^AIqIi'^ and circular velocities 
of Vcirc ~ 100 kms~^. The galaxy in the left panel ex- 
hibits a smaller scatter in the local TFR than the galaxy 
in the right one. By comparing the recent formation his- 
tories of these galaxies, we see that the large scatter in 
the second one can be associated to strong interactions 
(Pinter' ~ 0.4 and f;^™ - 0.2) taking place during the last 
6 Gyrs and which ended up in a recent important merger 
event (Fmerger ~ 0.25). Before this epoch, the K-ot/V^irc for 
the main progenitor remains close to ~ 1. In the case of 
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si.o/K:irc, the behavlour is similar but with more smooth 
variations. Note that simultaneously with the interactions 
and merger events, the gas fraction increases and the per- 
centage of gas in the disc component decreases indicating 
the presence of gas inflows which contribute more impor- 
tantly to the formation of a spheroidal component. In the 
case of the galaxy in the left panel, the lack of mergers 
and interactions during the last 5 Gyr is in agreement with 
the small variations of K-ot/K:irc and the high values of 
-^disc- Note, however, that its main progenitor experiences 
an increase of Vi-ot/Kirc 7 Gyrs ago, during a merger event 
(-^merger ^ 0.15), while 9 Gyr ago, a strong interaction 
(-^inter' ~ 0.6) Can be linked to a decrease by more than 0.1 
dex in K-ot/K;irc- Once again, the trends for si.o/K:irc are 
similar to those appreciated for Vrot/V'circ but with a more 
smooth evolution. Note that, as seen before, -Ffjgj. reaches 
a minimum during the period of strong interactions. 

Finally, in the lower panels, we compare the evolution 
of galaxies with low stellar masses of w 10^-'^MQh~^ 
and circular velocities of Vcm ~ 50 kms~^. The galaxy 
in the left panel shows a smaller scatter in the local TFR 
than the one in the right panel. Even in these small galactic 
systems, the scatter in Vrot/Vdrc seems to be related to the 
presence of merger events. The galaxy in the left panel has 
not suffered significant mergers or interactions during the 
last 5 Gyrs in agreement with its small changes of l^ot/K;irc- 
Note, however, that the progenitor in the main branch has 
been affected by mergers and interactions 8 Gyr ago, which 
led to a peak in the SFRe and caused changes by ~ 0.25 
dex in T4ot/K;irc and < 0.15 dex in the case of si.o/Kirc- For 
the galaxy in the right panel, the distortion of the rotation 
curve at z ~ can be again related with recent important 



mergers and interactions: F-^^^\ ~ 0.3 and 



merger 



0.3. 



Therefore, our results suggest that the hierarchical ag- 
gregation of the structure strongly influences the evolution 
of the galactic gas component within a A-CDM Universe, 
producing TFR outliers in concordance with previous 



works (e.g. Simard fc Pritchct 1998; Bar ton et all 1200 it 



KannaoDan fc Bartonii2004: iBohm et al.ii200^ Flores et al.l 



2006; Weine r et all 120061: iKassin et al.ll2007D . Mergers and 
interactions seem to be important at regulating the gas 
kinematics and star formation process for the whole range 
of stellar masses covered by these simulations, affecting 
the evolutionary tracks of systems on the TFR-plane. All 
these phenomena seem to be responsible for the scatter we 
foun d in the TFR. As p reviously reported in the literature 
(e.g. IKassin et all |2003), we obtained that the combina- 
tion of Vrot and a leads to a velocity scale more stable 
during merger and interaction events reducing the scat- 
ter in the sTFR and bTFR. We obtained, for example, 
that I log(K-ot/Veirc)| < 0.5 for the whole sample of sim- 
ulated systems while | log(si.o/Vcirc)| < 0.3. Note, how- 
ever, that some observational wor ks reported scatter s in 
the sTFR greater than ~ 1 dex (e.g. IKassin et al.l[2007l) . In 
this context, it is worth mentioning that the dispersion in 
the observed TFR relation might depend very sensitively 
on t he observational a pproach and modelling techniques 
(e.g. iMiller et al.ll201lD and it is yet not clear which is the 
intrinsic scatter of the rela t ion. B y analysing numerical sim- 
ulations, [Covingtoner^ (|2010t ) found scatters larger than 
1 dex. Nevertheless, those simulations were pre-prepared 
mergers of two isolated galaxies. In this work, by studying 



cosmological simulations, we are able to follow the evolu- 
tionary history of galaxies in a self-consistent way. 



4.2. Fingerprints of formation liistories on the TFR-plane 

In this section, we investigate how the assembly histories of 
galaxies in a cosmological scenario induced changes in the 
TFR-plane searching for fingerprints of particular events 
which, then, can be used to understand observations. 

In Fig. ini we can appreciate the evolution on the sTFR 
and bTFR planes of the six galaxies described in the previ- 
ous section. If the TFR is constructed by using V^irei it pro- 
vides us with information about how the stellar and bary- 
onic masses assembled within the potential well of the dark 
haloes hosting the galaxies, while if Vrot is used, it gives 
information related to the gas kinematics. To assess when 
galaxies are outliers of the mean TFR at different redshift, 
we also include the mean sTFR and bTFR of the who le pop- 
ulation as a function of redshift taken from de Ross i et all 
(201(|). For the sake of clarity, we did not include in Fig. 
|5| the tracks given by si.o but they can be easily derived 
from previous section: si.o would in general lead to similar 
trends to those obtained from V^ot but with a less noisier 
evolution and, given that si.o is larger than Kot by defi- 
nition, the tracks would be also displaced towards larger 
velocities. 

By comparing Fig. |S] and Fig. |21 we can understand the 
origin of the changes as the systems evolve. In particular, 
we see that mergers lead to an increase of Vdrc and mass 
in different amounts. With respect to interactions, in some 
cases, we detect a decrease in Vdrc due to interactions while 
in others the interaction generates an increase of the stellar 
and baryonic mass with an almost constant V^ire- When 
following the TFR tracks using Vrot we found a very noisy 
evolution which is affected considerably by interactions and 
mergers that produce outliers of the TFR. As expected, we 
see in Fig.jHjthat the tracks given by V^irc are smoother than 
those given by V^ot ■ Most mergers and interactions generate 
a scatter in the TFR greater than the level of evolution 
that we estimated since 2 = 3 for the relation based on 
Vciic in these simulations. As we noted before, these events 
can shift the T^ot towards lower or greater values than V^ire- 
These findings suggest that, from the observational point 
of view, selection effects and noise are likely to mask the 
actual evolution of the TFR. 

In the case of the galaxy corresponding to the upper left 
panel, the most important evolution in the TFR-plane took 
place more than 8 Gyr ago, during an epoch of important 



interactions and mergers (i^j 



2fivir 



4, F„ 



0.3) which 



generate an increase by more than 1.5 dex in stellar and 
baryonic mass and by ^ 0.4 dex in circular velocity. 

In the case of the galaxy in the upper right panel, the 
track given by Vrot approximates the one given by Vcirc 
with the only exception of the point associated to a look- 
back time of ~ 9 Gyr when the system experienced an 
important merger event (Fmerger ^ 0.25). As a consequence 
of the merger, Vrot decreases by ^ 0.4 dex departing from 
Kire and generating a TFR outlier. After ^ 1 Gyr, the disc 
component was reconstructed, with Vrot approaching V^ire- 

With respect to the galaxies in the middle panels, we see 
that the evolutionary tracks given by VJ-ot are noisier than 
those of more massive galaxies. This suggests that inter- 
actions and mergers in these simulations are more efficient 
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at producing instabilities in the disc component of systems 
with shallower potential wells. In these galaxies, the varia- 
tions in the rotation velocity between two time steps of the 
simulation are, in some cases, larger that the total level of 
evolution of the TFR based in Vdrc for 3 < z < 0. In par- 
ticular, the galaxy in the middle right panel reaches z ^ 
as an outlier of the TFR as a consequence of recent impor- 
tant interactions followed by a merger event (F^J^^^' ^ 0.4, 

-'^inter ^ O'^i -Plnorgcr ^ 0.25). In this figure, it is evident 
that during mergers and interactions the tracks given by 
Vcirc and Vrot may diverge. 

In the lower panels of Fig.[9j we can appreciate that low- 
mass galaxies also describe a similar behaviour. The galaxy 
in the left panel becomes an outlier of the TFR given by 
Vrot at a lookback time of ~ 9 Gyr due to interactions and 
mergers. After ^ 1 Gyr, this galaxy recovered its rotational 
equilibrium with VJ-ot approaching Vcirc again. In the case of 
the young galaxy in the right panel, the presence of mergers 
and interactions displaced this galaxy away from the mean 
TFR based in Vcirc- We can also appreciate that due to a 
recent merger event, Vcirc and Vrot evolve in opposite ways. 

Although we have pointed out only the correlation be- 
tween mergers and interactions, from the previous section, 
we know that outflows and inflows can also be associated to 
sharp changes on the TFR-plane. However, in most cases 
they also occurred close to mergers or interactions as a 
second order effect since these violent events have been 
proven to be efficient at triggering g as inflows and star- 
bursts (which can drive outflows, e.g. Barnes fc Hernauisj 
[ToollMihos fc HernauistilT996tlTisserall2000D . 



4.3. Statistical analysis 

Our results for these six galaxies suggest that the hierar- 
chical building up of the structure strongly influences the 
evolutionary paths of galaxies in the TFR-plane modulat- 
ing the evolution of its scatter with cosmic time. In or- 
der to generalise these findings, we performed a statisti- 
cal study by extending the previous analysis to the whole 
sample of simulated galaxies. Within the hierarchical ag- 
gregation paradigm, galaxies are afi^ected by different kind 
of physical processes such as mergers, interactions, star for- 
mation, outflows, gas accretions, etc. All these processes act 
together to drive galaxy evolution and hence it is not easy 
to disentangle the effect of each individual event on the vari- 
ation of the properties of galaxies. It is also possible that 
the combined actions of two or more different astrophysical 
processes compensate each other to generate a insignificant 
change of the mean properties of the system. Therefore, we 
decided to focus our attention on galaxies which have expe- 
rienced extreme physical events in the very recent past (in 
the previous 2 Gyr of galaxies selected at a given z). For 
these cases, we expect the evolution of the mean properties 
of galaxies to be mainly affected by a particular extreme 
physical event. 

We assumed that a galaxy has experienced an ex- 
treme event when one of the following conditions is true: 

-^m erger 

> 0.3 (Fig. [101 upper left panel), F.^^^[ > 0.6 (Fig. 
[TUl upper right panel), AMbar/Mbar < -0.3 (Fig.[Tni lower 
left panel) and A/gas > 0.3 (Fig. [TUl lower right panel). 
AMbar/-Mbar IS defined as the change in the baryonic mass 
of the galaxy between two time steps of the simulation if the 
time interval is not greater than 2 Gyr and is normalised to 



the initial baryonic mass. In these simulations, a decrease 
in the baryonic mass of a galaxy (AMbar/Afbar < 0) im- 
plies a significant loss of gas and, therefore, the condition 
AMbar/Mbar < —0.3 allows US to select systems which have 
been subject to strong outfiows events. In fact, these simu- 
lated galaxies would have lost at least 30% of their baryonic 
mass during the last 2 Gyr. A/gas is defined as the total 
variation in the gas fraction of the galaxy with time. In 
these simulations, the gas inside the galaxy tends to de- 
crease because of the star formation process or ejections by 
galactic winds. Hence, the gas mass can only increase by in- 
fall of surrounding material and/or by mergers. Therefore, 
our condition A/gas > 0.3 allows us to identify galactic 
systems which have suffered important gas accretion. 

We estimated the distribution of the absolute variations 
of log(K-ot/K:irc) and log(si.o/V^irc) of galaxics which have 
suffered extreme events, as shown in Fig. 1101 For compar- 
ison, we calculated the corresponding distribution for the 
whole sample of simulated galaxies. From the histograms, 
it is clear that the distribution of galaxies which have 
been subject to extreme events are biased to larger val- 
ues of |Alog(V[.ot/K;irc)| than those of the whole galaxy 
sample, indicating that all these processes tend to gen- 
erate TFR outliers by disturbing the gas kinematics. In 
particular, most events (from any of the mechanisms con- 
sidered) induce mean deviations of ~ 0.1 dex, with maxi- 
mum offsets in the range [0.4,0.5] dex. We can also see that 
|Alog(K-ot/Kirc)| < 0.5 but, if we use si.o instead of K-ot, 
the distribution is narrower with |Alog(si.o/T4irc)| < 0.3, 
approximately. In the case of |A log(si.o/T4irc)|, disturbed 
galaxies tend to be biased towards larger values when asso- 
ciated to mergers and interactions. However, the deviations 
of s 1.0 from Vdic are less significant in the case of important 
infiows and outfiows. This suggest that si.o can account 
better for dispersion associated to these later effects than 
to those which could be imprinted by more violent effects 
such as mergers and interactions. 

In order to analyse if there is any trend for some of the 
mechanisms to drive a positive or negative variation of the 
rotation velocity relative to the circular one, we estimate 
the percentage of galaxies with A log(T^ot /Khc) < in each 
of the subsamples. We found a weak trend for galaxies ex- 
periencing important mergers or outflows to have negative 
variations (61% and 58%, respectively) while strong inter- 
actions and gas inflows can drive either positive or negative 
outliers on the TFR-plane (49% and 53%, respectively). 

These results indicate that the joint action of merg- 
ers, interactions, starbursts, outflows and gas accretion can 
affect the TFR producing variations in log(K-ot/Vi3irc) as 
larger as ~ 0.5 dex, at least in these simulations. Hence, 
at a given stellar or baryonic mass, these events can pro- 
duce a scatter in the TFR-plane larger than the mean level 
of velocity evolution since z ^ 3 (^ 0.1 dex). In fact, for 
any of the subsamples analysed in Fig. [TOl while ^ 25% of 
the galaxies have variations in log(14ot/Kirc) of ^ 0.1 dex, 
55% of the systems exhibit values > 0.1 dex, which might 
mask the evolution of the TFR in observational studies. 
However, as we have seen, by combining V^ot and a in the 
deflnition of the kinematical indicator, it is possible to gen- 
erate a velocity scale more robust against disturbances of 
the gas kinematics being able to reduce the scatter in the 
simulated sTFR and bTFR by a factor of - 2. 
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5. Conclusions 

We studied the evolution of the gas kinematics of galaxies 
by using numerical simulations which include a physically 
motivated SN feedback. We focused on the origin of the 
scatter of galaxies in the TFR-plane and the connection 
with the formation histories of these systems. In particu- 
lar, we analysed the role of mergers and interactions on 
the determination of gas kinematics and of the features 
of rotation curves during the assembly o f galaxies. We ex- 
tended the work bv lde Rossi et al.l ()2010D . who studied gas 
rotation-dominated systems, and explore the whole sample 
of galactic systems including also gas dispersion-dominated 
galaxies. It is worth mentioning that our simulated sam- 
ple does not contain realistic spiral-like galaxies as, in gen- 
eral, simulated systems have dominant stellar spheroids. 
Nevertheless, all analysed galaxies have surviving gaseous 
discs which are found to be good tracers of the potential 
wells. 

Our main results can be summarised as follows: 

— Our simulations support the claim that V^ot is the best 
proxy for Vdrc at i?max for galaxies of all morphol- 
ogy types. For gas rotation-dominated systems, V^ot is 
a good proxy of Vdrc at all radii. As the dispersion- 
dominated component increases, cr/Vrot gets larger, spe- 
cially in the outer parts of the systems. This leads to an 
increase of the scatter of the TFR when considering sys- 
tems with more disordered gas kinematics. Nevertheless, 
for systems with large dispersion-dominated compo- 
nents, Vrot is still a good proxy for Vdrc at i?max- In par- 
ticular, for spheroidal dominated galaxies, we obtained 
a mean cr/Kot of ~ 0.7 at J?max in good ag reement with 



recent observations ([Catinella et al.ir2012 ). 
The use of the kinematical indicators sp. 5, si.p and S 
(IWeiner et al.ll2006HKassin et al.l[2007t ICovington et al.l 
in our galaxy sample, leads to a reduction of the 
scatter of the TFR in all cases. Moreover, all of them, 
together with Vrot , lead to the tightest relations if eval- 
uated at i?max- In particular, by comparing so.5 and 
S with K;irc, we obtained that the former underesti- 
mates K;irc, while the latter overproduces it. Although 
in the literature sq.s has been used more frequently, we 
find that, in the case of these simulations, si.o is a bet- 
ter kinematical indicator as it does not only reduce the 
scatter of the TFR but it is also a good proxy of Vch-c 
at large radii (r > 0.5-Rbar)- By restricting our study 
to subsamples with different baryonic D/T or differ- 
ent relative angular momentum orientation between the 
gas and stellar phases, we obtained, on average, similar 
global trends. Nevertheless, galaxies with misalignments 
between the gaseous and the stellar angular momenta or 
with lower baryonic D/T present larger scatter around 
the mean behaviour. 

Mergers and interactions within a A— CDM scenario 
can regulate the star formation process as well as drive 
inflows and outflows of gas affecting the evolutionary 
tracks of galaxies on the TFR-plane and generating 
TFR outliers. In particular, we obtained that merger- 
induced velocity disturbances in our cosmological simu- 
lations are typically smaller (< 0.5 dex) than predicted 
by pre-pre pared mer ger simulations of Milky- Way type 
galaxies (|Covington et al.. 2010). This is a consequence 
of the fact that the cosmological merger histories of the 
galaxies in the simulated mass range involved typically 



more minor mergers than the major mergers commonly 
built in merger simulations. 

— We detected that mergers and interactions can generate 
either an increase or decrease of K-ot/14irc- Nevertheless, 
statistically, we obtained a weak trend for galaxies sub- 
ject to important mergers to be biased to negative vari- 
ations (~ 61%), while strong interactions can generate 
either positive or negative changes with similar prob- 
ability. We also found that gas infall or outflows can 
lead to TFR outliers. In particular, the statistical study 
shows a weak trend for outflows to produce negative 
variations 58%) while gas inflows can drive either 
positive or negative changes. 

— By comparing the distributions of the absolute vari- 
ations of log(si.o/V;ii.c) and log(yrot/K:irc) during the 
evolution of all simulated galaxies, we obtained that, for 
the former case, the scatter is reduced by a factor of ^ 2, 
exhibiting a narrower distribution (|Alog(si.o/K;irc)| < 
0.3). In particular, mergers and interactions can be as- 
sociated to a larger scatter in the TFR even when using 
si.o; conversely, the effects of inflows and outflows seem 
to be more efflciently accounted by this kinematical in- 
dicator. 

— According to our results, extreme physical events in sim- 
ulated galaxies (e.g. mergers, interactions, outflows and 
inflows) lead to a scatter on the TFR-plane larger than 
the mean level of evolution of the TFR since 2; ~ 3. In 
addition, by enhancing the star formation rate, merg- 
ers and interactions also affect the luminosity and M/L 
ratio of the systems, which can induce observational off- 
sets from the TFR. Therefore, much work is still needed 
from the observational and theoretical point of views to 
address these problems. 
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Fig. 8. Formation histories for six typical simulated galaxies at 2: = in the simulations, with different stellar (M*) 
and baryonic (Mb) masses and circular velocities (Kirc)- The masses are given in logarithm and in units of Mq/i"-^ The 
circular velocities are given in logarithm and in units of kms~^. In panels A, we show the evolution of F-^^^l (dotted 

line) and F^J^^^' (dashed line). In panels B, we can appreciate Fmoigcr- Panels C show the time evolution of the gas 
fraction (dotted line), baryonic mass relative to the final one (triple-dot-dashed line) and F^^^^ (dashed line). Panels D 
show the star formation efficiency SFRe in yr~^ and, in panels E, we can appreciate the ratio V^ot/Kirc (dashed lines) 
and si.o/K;irc (solid lines) evaluated at i?bar- See the text for definitions of variables. 
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Fig. 9. Evolutionary tracks of the six galaxies shown in Fig. |8] on the sTFR and bTFR planes. The background grey 
lines depict the mean sTFR and bTFR obtained by using Vi^irc at i?bar as a kinematical indicator at z = 3,2,1,0.4,0 
(dotted, dashed, dot-dashed, triple-dot-dashed, solid line, respectively). The black lines indicate the track of the given 
galaxy on the TFR plane by using V^irc (solid line) and V,ot (dashed line) at i?bar as the kinematical indicators. In each 
panel, there is an arrow indicating in which direction the galaxy evolves along each path. 
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Fig. 10. Distribution of the absolute variation of log(K-ot/Kirc) (main panels) and log(si.o/Kirc) (insets) during the 
last 2 Gyr of evolution for galaxies which have been subject to extreme processes (dashed lines): Fmcrgcr > 0.3 (upper 
left panel), F-^^^',. > 0.6 (upper right panel), AMbarZ-^bar < —0.3 (lower left panel) and A/gas > 0.3 (lower right panel). 
Solid lines show the corresponding distributions for the whole sample of simulated galaxies at available redshifts. 
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